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We study nature of singularities in anisotropic string-inspired cosmological models in the presence 
of a Gauss-Bonnet term. We analyze two string gravity models- dilaton-driven and modulus-driven 
cases- in the Bianchi type-I background without an axion field. In both scenarios singularities can 
be classified in two ways- the determinant singularity where the main determinant of the system 
vanishes and the ordinary singularity where at least one of the anisotropic expansion rates of the 
Universe diverges. In the dilaton case, either of these singularities inevitably appears during the 
evolution of the system. In the modulus case, nonsingular cosmological solutions exist both in 
asymptotic past and future with determinant D = +oo and D = 2, respectively. In both scenar- 
ios nonsingular trajectories in either future or past typically meet the determinant singularity in 
past /future when the solutions are singular, apart from the exceptional case where the sign of the 
\ time-derivative of dilaton is negative. This implies that the determinant singularity may play a 

crucial role to lead to singular solutions in an anisotropic background. 
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>* ' I. INTRODUCTION 

Superstring theory continues to be of interest as a possible candidate to unify all fundamental interactions in nature 
Jy. It is known that there exist five supersymmetric perturbative string theories which are classified as the type I, type 
IIA, type IIB, SO (32) heterotic and E$ x _E g heterotic string. Recently it was found that these theories are connected 



O 
O 



by dual symmetries, which leads to the conjecture that each theory appears as one of five branches of a unified 



theory, called M-theory In particular Hofava and Witten || showed that the 10-dimensional E$ x E$ heterotic 
^Jy string theory is equivalent to an 11-dimensional M-theory compactified to M 10 x S 1 /^. Then the 10-dimensional 
L* ' spacetime is expected to be compactified into M 4 x CY 6 , in which case the standard model particles are confined 
on the 3-dimensional brane. This gives rise to the well-known brane world scenario [Q where the extra dimension is 
5^ ■ noncompact and gravity is effectively 3-dimensional. 



For cosmologists it is very important to test the viability of string theories by extracting various cosmological 
implications from them One of such attempts is the pre-big-bang (PBB) scenario || based on the low energy 
effective action of string theory. In this scenario there exist two branches of solutions by assuming a T-duality. One 
of which (t < 0) corresponds to the stage of pole-like inflation driven by the kinetic term of the dilaton field, and 
another (t > 0) is the stage where the curvature continues to decrease. However it is difficult to smoothly connect 
these two branches without singularity in tree-level string action jl],|8| . 

One is required to take into account quantum loop or derivative corrections in order to overcome such singularity 
problems. In fact, Antoniadis, Rizos, and Tamvakis Q included a Gauss-Bonnet term to the tree-level string effective 
action with dilaton and modulus fields, and showed the existence of nonsingular cosmological solutions. In this case 
nonsingular behavior of solutions is mainly determined by the evolution of the modulus field. Therefore allowed ranges 
of parameters were analyzed in the absence of the dilaton in the flat Friedmann- Robertson- Walker (FRW) background 
fiof (sec also ref. Jll|). Since it is important to confirm the generality of singularity avoidance even starting from an 
anisotropic spacetime, several authors analyzed nonsingular cosmological solutions in the Bianchi type-I spacetime 
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without dilaton (lJJT^l and with dilaton |TJ] . The presence of the modulus coupled to the Gauss-Bonnet term allows 
the existence of nonsingular solutions unless the dilaton controls the dynamics of the system. 

In order to understand how nonsingular or singular solutions appear, it is necessary to classify nature of singularities 
in an anisotropic background. In particular, the main determinant D of the system is an important quantity to describe 
the singularities. When only the dilaton field <f> is coupled to Gauss-Bonnet term in the Bianchi I background, it was 



conjectured in ref. |15| that nonsingular cosmological solutions in future crosses the determinant singularity (D = 0) 
in past when cj) is positive. While this singularity was found more than ten years ago Jl6[ , only now we begin to 
understand its importance. The similar kind of singularity also appears in the context of black hole inner solutions 



in the presence of dilaton coupled to gravity via the Gauss-Bonnet term 17 lqj. In this paper we shall make detailed 
analysis about nature of singularities both in dilaton- and modulus-driven cosmologies in the Bianchi I background. 
We do not include an axion field in our analysis, but it is important to emphasize that its effect is generally vital as 
studied in Refs. JlSj|. We will classify other kinds of singularities where at least one expansion rate diverges. These 
investigations are important to understand how nonsingular solutions emerge in the modulus-driven case. In addition 
our analysis will be useful to construct more complicated nonsingular string-inspired models in the presence of other 
fields such as axion. 

This paper is organized as follows. In Sec. II we show background equations in anisotropic string-inspired models 
with dilaton or modulus fields. In Sec. Ill we study nature of singularities in dilaton-driven cosmology both for 
positive and negative <fi cases. Sec. IV is devoted to the modulus-driven cosmology where both of nonsingular and 
singular solutions exist. We present summary and discussions in the final section. 

II. THE MODEL AND BACKGROUND EQUATIONS 

We begin with the action |j] |l4| . 



S = Id x^J—g 



ii?-i(V0) 2 + /(0)i^ B 



(2.1) 



written in the Einstein frame. Here R is the scalar curvature and <f> denotes a scalar field which is either dilaton or 
modulus. /(</>) depends on string theories, whose explicit forms are given later. We do not consider the multi-field 
system of dilaton and modulus fields jl4j induced from the one-loop effective action of heterotic string theory. In 
addition we neglect the anti-symmetric tensor H^ v \ and the curvature terms higher than the second order. The 
Gauss-Bonnet term, i?GB, is defined as 

R 2 GB = R 2 - AR^R^ + R^ af) R^ af 3. (2.2) 

We normalize time and spatial coordinates by the string length scale ^/X^ as x^ — x^/^/X^, and the scalar fields as 
(\> = </>\/A7. Hereafter we drop bars for simplicity. 

Let us consider the Bianchi type-I spacetime whose metric is given by 

ds 2 = -dt 2 + a 2 (t)dx 2 + b 2 (t)dy 2 + c 2 (t)dz 2 , (2.3) 

where a(t) , b(t) , c{t) are the scale factors in an anisotropic background. We define the anisotropic expansion rates, 
p(t),q(t),r(t), as 

P{t) = -, q(t) = b r, r(t) = - (2.4) 
a b c 



2 



where a dot denotes the derivative with respect to t. It is also convenient to introduce new variables h,a,0: 

p(t)=h + a + V3(3, q(t) = h + a-V3/3, r{t) = h-2a. (2.5) 

Here h = (p+q + r)/3 is an average expansion rate, which is the generalization of the Hubble parameter in an isotropic 
case, and a and [3 correspond to parameters of anisotropy. The triangle in (a, /3)-plane, 



a + V3/3>-h 7 a-V3(3>~h 7 a < h/2 , 



(2.6) 



which extends around the isotropic point a = /3 = 0, represents the regions where the Universe expands in all 
directions. In the outside region the Universe is contracting at least in one direction. 
The dynamical equations for the background are written as 



(1- 


f 8rf)(q 


fi 2 ) 


+ (l + 8«/)(r + 


r 2 ) + (l + 8/)gr + 


2 V 


-0, 


(2.7) 


(1- 


f8r/)(p 




+ (l + 8p/)(r + 


r 2 ) + (l + 8/>p4 


2 V 


= 0, 


(2.8) 


(1- 


f8g/)(p 




+ (l + 8p/)(g + 


g 2 ) + (l + 8/)pq4 


2 V 


-0, 


(2.9) 




- (p + q + 




— 8/' {pgr + pgr 


+ pgf + pqr(p + q 


fr)} 


= 0, 


(2.10) 



together with the constraint equation, 

1 •„ 

pq + qr + rp + 2Apqrf — —(f) = . 
Defining a 4-dimensional vector, x = (p, q, r, ifi), eqs. (2.7)-( 2~To| ) can be written in the matrix form, 

= y , 



(2.11) 



(2.12) 



where y = y(p, q, r, < 



and 



Z = 



l + 8/>r 1 + 8/^3 8/V 

l + 8f<j>r l + 8/>p Sfrp 

l + 8f'H l + 8/'# 8f P q 

-Sf'qr -Sfrp -Sf'pq 1 



(2.13) 



The determinant of Z yields 

D = 2 + 16f<Xp + q + r) - QAf' 2 {p 2 q 2 + q 2 r 2 + r 2 p 2 ) + 128/' 2 pqr(p + q + r) + 128 f' 2 <j> 2 {pq + qr + rp) 

+ 1024 j '' 3 4>pqr(pq + qr + rp + <j) 2 ) + 12288 j ~' 4 4> 2 p 2 q 2 r 2 . (2.14) 



In the case of D ^ 0, the solutions of eqs. (2.7)-(2.10) are given by x = Z~ x y. When D vanishes, however, we cannot 
proceed numerical calculations further. This "determinant singularity" plays an important role in the anisotropic 
background p5[. 



From Eq. (2.11) we find the constraint 



'^2 



in which case d> is solved as 



a 2 + f3 2 < h 2 + m{pqrf) 



24pqrf ± y/(2Apqrf') 2 + 2(pq + qr + rp) 



(2.15) 



(2.16) 



When / — anisotropy parameters are restricted in the circle, a 2 + 1 < h 2 . If the Gauss-Bonnet term is taken into 
account, we have the wider allowed range of anisotropy parameters given by eq. ( 2.1 5| ) . 
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III. DILATON-DRIVEN CASE 



Firstly we consider the dilaton-drivcn case with 



-2rf> 



(3.1) 



where the string coupling, A, takes a positive value. We set A = 1 in our numerical simulations. In the scenario (3.1) 
nature of singularities was analyzed in rcf. EH] in the Bianchi I background in the case of the plus sign in the rhs of 
Eq. (2.16). Hereafter we shall make detailed analysis about the property of singularities in both signs of eq. f2~ll| ). 

The asymptotic behavior of solutions in past and future can be analyzed by assuming the following power-law forms 
for the expansion rates: 



p = Cl \t\ s , q = c 2 \t\ 



oai*r 



(3.2) 



In order for the / term in eqs. (2.7)-( 2.10 ) to have a power-law dependence, the dilaton is required to take the form 



4> = 4>o + c 4 ln|i| 



(3.3) 



When the contribution from the Gauss-Bonnet term is negligible, one has (f> 2 — 2(ciC2 + C2C3 + C3C1 ) |i| 2s from 



eq. ( 2.11 ). Comparing this with eqs. (2.7)-( 2~Tc| ) and we find 

s = -l, ci + c 2 + c 3 = sign(i), 



2,2,2 
Ci + c 2 + c 3 



c2 = 1, 



C4 = 2(cic 2 + c 2 c 3 + C3C1) , 



(3.4) 



In the absence of the dilaton (C4 = 0), the solution ([O]) for t > represents the vacuum Kasner solution where 
the Universe is expanding in two directions and contracting in one direction. The interpretation of this solution is 
that large anisotropies are required at least in one dimension as t — > 0, in order to make the spacetime curved by 
anisotropies. 

The situation is different when the dilaton is taken into account. For example, when \/l/2 < c 4 < \/2/3 the 
Universe is expanding in all directions for t — > 00, while there exist both Friedmann- and Kasner- type solutions 
for < c 4 < •v/1/2. Note that in the limit where the Gauss-Bonnet term is negligible (/' — > 0) the determinant 



approaches a constant value D = 2 from eq. (2.14) 



When the Gauss-Bonnet term is dominant in eq. (2.11) (\pq + qr + rp\ <§; |24p«jr/|), one has 



-6Xcic 2 c 3 \t\ 6s e 



3s 20 



by using (3.2). Integrating this equation with respect to t, we easily find that 4> = sign(t) (3s + l)/(2|f|). Combining 



this with eqs. ( 2.7 )-( 2.1C ) gives 



-2, 



c 4 



cic 2 c 3 = sign(i) 



5e 20o 

12A 



(3.5) 



Since C1C2C3 < for t < 0, the Universe is either contracting in all directions or expanding in two directions and 
contracting in one direction. When t > 0, the Universe is either expanding in all directions or expanding in one 
direction and contracting in two directions. For the asymptotic solution (|3.5| ), we have that <ft oc and /' oc |t| 5 , 



for \t\ 



In spite of this divergent behavior of the 



for |t| -> 00. 



in which case the determinant is given by D oc |i| 6 — * 00 
determinant, the solutions are nonsingular with p,q,r oc |i| 

We shall classify the cases where the solutions of eqs. ([2~7|)-( 2.1C ) exhibit singular behavior. When the system 
passes through the determinant singularity (D = 0), eq. ( |2.12j ) indicates that x = (p,q,r,(f>) diverge. This singularity 
appears in an anisotropic background where three expansion rates are multiple- valued functions of time |^3| . It is 
also a physical singurality where the curvature invariant, Rv a P R^ap, diverges due to the divergence of the time 
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derivative of the expansion rates. Near the determinant singularity, the expansion rates and the scalar field can be 
expanded as [n5| 



hi = h is + hny/\t - t s \ + h l2 (y/\t - t s \) 2 + 



= S + MV\t-ts\) 2 + MV\t-ts\) 3 + ■■■, 



(3.6) 
(3.7) 



where hi — p,q,r (i = 1,2,3), and t s is the time at singularity. This means that hi and <j) diverge as t — > t s , while 
p, q, r, <f> are finite. This property is different from the ordinary kind of singularity where p, q, r do not stay finite. The 
determinant singularity plays a crucial role in an anisotropic background. 
The ordinary kind of singularities can be classified as 



(i) p~p /(t-t s ), q 

(ii) p~p Q /(t-t s ), q 

(iii) p~Po/(t~t s ), q 



q , r ~ r , 
p /(t-t s ), r~r , 
Po/(t-t s ), r ~ Po/(t-t s ) , 



(3.8) 
(3.9) 
(3.10) 



-0) one has p 



-oo, 



where po,qo> an d r$ are constants with po > 0. If the time-direction is futurewards (t — t s 
while p — > oo for t — t s — ► +0. The asymptotic forms of the determinant (2.14) depend upon the cases presented 
above. 



In the case (i) with a plus sign in rhs of Eq. ( 2.16| ) , the signs of qo and ro are the same and the asymptotic form of 
tj) is given by <f> ~ (q + r ) / (3Xe~ 2 ^ s q r ) for t — t s — * +0, with (j> s being a constant. Then the determinant D yields 



from eq. (2.14): 



D 



"Pi 



ro 



oo 



(for q Q r > 0) . 



In the minus sign of eq. ( 2.16| ), one has D — > — oo for go r o < and t — t s 



In the case (ii) with a plus sign of of eq. ( 2.16 ), one has </> ~ l/(3Aroe 2 ^ s 
determinant should asymptotically takes the form 



-0. 

) for ro > and t — t s 



(3.11) 



-0. Then the 



D 



oo 



(for r > 0) . 



(3.12) 



3(<-i s ) 4 

In numerical analysis we did not find this case for the parameter ranges and initial conditions we adopt. Typically 
ordinary singularities are dominated by the cases where one or three expansion rates tend to diverge. When one 
chooses the minus sign of eq. ( 2.16 ), the asymptotic behavior is D — > — oo for r Q < and t — t s 
In the case (iii) with a plus sign of of eq. ( 2.16|) , we find <j> ~ t/(Xp e~ 2 ^ 3 ) and 



-0. 



D 



0, 



(3.13) 



for t — t s — > +0. Strictly speaking this holds only for the isotropic case (a = /3 = 0) where all expansion rates are 
the same. In this case the solutions do not cross the determinant singularity, although they approach D = as 
t — t s — * +0. When small anisotropies are included, the trajectories can pass through D = 0. For the minus sign of 



of eq. ( 2.16 ) the asymptotic behavior of the determinant is not described by eq. ( 3.13 ), as we will see later 
In another limit, t — t. 



—0, the signs of diverging expansion rates in eqs. ( |3.8| )-(3.10) are reversed, in which case 
the asymptotic forms of D are altered. Nevertheless the determinant of the cases (i) and (ii) generally approaches the 
asymptotic value, D = — oo. We shall confirm this by numerical investigations in subsequent sections. 
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A. Plus sign of eq. (2.16) 



We first analyze the case of the plus sign in eq. ( 2.16| ). When is largely positive, the term /' = — (A/8)e 



is negligible in (2.16), implying that <j) is positive as long as pq + qr + rp > 0. In this case <p increases toward the 
future, which results in nonsingular asymptotic solutions with determinant, D rs 2. When we go back to the past, 
there are two possibilities for the evolution of the determinant. One is the case where the solutions pass through the 
determinant singularity (D — 0) and another is the one where D goes toward infinity (D — > +oo). Our numerical 
investigations suggest that the latter case does not occur for the plus sign in eq. ( 2.16 ), which means that solutions 
nonsingular in future (D — > 2) meet the determinant singularity in past, irrespective of the initial values of cj). 

In the left panel of Fig. |l| (a) we plot the evolution of the expansion rates for the anisotropy parameters a = 
0.05, (3 = 0.05 with initial conditions h = 0.16 and (j> — 0, corresponding to D w 1.7 at t = 0. The determinant 
continues to grow until it approaches the finite value D = 2 as t — > oo, which indicates that this solution belongs to 
the nonsingular solution in future given by eq. (3.4). Note that for t > the Universe is expanding in all directions. 
When we solve the equations of motion pastwards (t < 0), the solution meets the determinant singularity around 
t = —0.82, thereby leading to the divergence of p, g,r and <f>. If we introduce a new time parameter, r, defined by 

r-/§, (3.14, 

it becomes possible to enter the region of the negative sign of D by overpassing the determinant singularity |l3| . This 
does not mean that we can remove the singularity by coordinate transformations. The determinant singularity is 
a physical one where the divergence of the curvature invariant is unavoidable even in other coordinates. Then the 
solution turns back futurewards and the determinant begins to decrease rapidly toward D — > — oo [see the right panel 
of Fig. [| (a)]. From Fig. [l] (a) we find that this belongs to the class of the case (i) with p ~ p , q ~ q ,r ~ r Q /(t — t s ) 
and poqo > [see eq. ( |3.1l| )]. The Universe is rapidly contracting in one direction (r — > — oo) for t — t s — > —0. As 
claimed in ref. Jig] , these trajectories can be understood as the pair creation of two branches (D > and D < 0) at the 
determinant singularity. For the nonsingular solutions in future, it is inevitable to cross the determinant singularity 
in past in an anisotropic background. Notice that in the isotropic case D is always positive and decreases toward zero 
as t — > — oo. In the presence of small anisotropies, however, the solution reaches the determinant singularity at finite 
past and fall into the D = — oo singularity as shown in Fig. |l| (a). 

When <p is large, the allowed anisotropy parameters in (a, (3) -plane lie inside a circle, given by 

a 2 +[3 2 <h 2 , (3.15) 

which comes from the constraint ( 2.15| ). In Fig. (a) we show nature of singularities in past and future in the (a, (3)- 
plane for h — 0.16 and (f> = at t = 0. In this case the allowed region is approximately described by eq. fl3.15p . We 
find that all solutions inside this region exhibit the determinant singularity in past (region II in Fig. ^) while they are 
nonsingular (D — > 2) in future (region I) . 

With the decrease of </>, the allowed region gets larger as found by eq. ( 2.15 ). Nature of singularities becomes more 
complicated due to the appearance of ordinary singularities. If the system is close to the isotropic case (|a|, \(3\ <C 1), 
one has D — > 2 for t — > oo and D — > at finite past as in the case of Fig. ^ (a). However, larger anisotropies alter this 
picture. For example, we show the evolution of the expansion rates in the left panel of Fig. |l| (b) for a — 0.15 and 
(3 = 0.05 with h — 0.16 and = — 1 at t = 0. In this case p is singular at finite past while q, r are positive constants, 
which means that D — > — oo from eq. (3.11). The solution exhibits the same kind of singularity (D — > — oo) in future. 
In this case we have numerically found the determinant is always negative (D = —1.02 at t = 0) and the solution 
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never crosses D = in both past and future [see the right panel of Fig. [| (b)]. In Fig. || (b) the asymptotic property 
of singularities is presented for h = 0.16 and <fi = — 1 at t = 0. There exists the region III where the determinant is 
singular (D — > — oo) in both past and future. We also find that the region IV with D — > +oo appears in past when 
anisotropy parameters are large. This corresponds to asymptotic past nonsingular solutions given by eq. (3.5) where 
quadratic curvature corrections are dominant. It is important to separate this case from the ordinary singularity with 
D = — oo, although this classification was not done in ref. [^5|. From Fig. ^ (b) the solutions with D = in future 
correspond to, in past, the determinant singularity or the regular solutions with D — > +oo. Notice that the parameter 
range of future nonsingular solutions are smaller compared to the case of Fig. |^ (a) . 

In Fig. U (c) we show the density plot for h = 0.16 and 4> — — 2 at t = 0. This suggests that future nonsingular 
trajectories are restricted to be very narrow near the isotropic point, a — = 0. We also find that the determinant 
singularity in future corresponds to the regular solutions with D — > +oo in past. As one example we plot in Fig. |l| 
(c) the evolution of p,q,r for h = 0.16 and <f> — —2 at t = with anisotropy parameters, a = 0.1, and (3 = 0.1. In 
this case the past asymptotic solution is categorized in the nonsingular solution given by eq. (3.5). Fig. || (c) shows 
that two expansion rates p, q are positive while r is negative for t < 0, implying C1C2C3 < in eq. (3J3). The solution 
comes regularly from the asymptotic past with D being decreased toward the future [see the right panel of Fig. [l] (c)]. 
It crosses the determinant singularity around t = 0.3, after which the determinant continues to decrease until the 
solution falls into the ordinary singularity with D — > —00. Note that p diverges as t — t s — > +0, while q and r approach 
positive constant values, in which case one has D — > —00 by eq. (3.11). This trajectory can be regarded as the pair 
annihilation of two branches with D > and D < 0. From Fig. || (b) and (c), the solutions nonsingular in past 
(D — > +00) meet the determinant singularity in future. In this case the solutions do not go futurewards beyond the 
determinant singularity. In Fig. |^ (c) the parameter regions with ordinary singularity in future (D = —00) typically 
correspond to those with the same singularity in past, in which case the asymptotic behavior of the expansion rates 
are the similar as in Fig. [I] (b) . Notice that in Fig. || (c) there exist some parameter ranges where the solutions meet 
the determinant singularity in both past and future. Although this case is very rare, it is still possible to cross the 
determinant singularity twice. 

From Fig. |we find that nonsingular solutions in future are of the determinant-type in past. This can be understood 
that the determinant evolves from D = 2 to D — 0, if we solve the equations of motions from asymptotic future to 
past [see the right panel of Fig. (a)]. We also analyzed other cases varying the values of h and (j> for 10~ 5 < h < 1 
and —15 < <f> < 15. For a fixed value of h there exists a minimal value of <fi which leads to nonsingular solutions in 
future [see Fig. |]. The allowed range of the dilaton for future nonsingular solutions gets wider with the decrease of 
h as shown in Fig. |§|. In all cases analyzed in our numerical simulations, the past singularity for future nonsingular 
solutions corresponds to the determinant-type in an anisotropic background. We also found that trajectories with 
nonsingular past asymptotic meet the determinant singularity in future. In addition we showed that the ordinary 
singularity with (3.8) appears, in which case the determinant is divergent as D — > —00. This is different from the 
case D — > +00 where the solutions are nonsingular in asymptotic past. In the next subsection we will analyze how 
the behavior of the determinant is altered for the different sign in eq. ( 2.16] ) . 
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FIG. 1: The evolution of p, q, r, (j>, and D in the dilaton-driven case for the plus sign in eq. ( 2.16 ). 
Each figure corresponds to initial conditions with h = 0.16 and (a) (j> = 0, a = 0.05, (5 = 0.05, 
(b) cj> — ~ lj Q = 0.15, P — 0.05, (c) (p = —2, a = 0.1, ft = 0.1, respectively. The determinant 
singularity (D = 0) can be passed through by introducing a new time parameter, r = J dt/D. 
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FIG. 2: Nature of singularities in (a,/3)-plane with initial conditions, (a) <f> — 0, (b) <f> = — 1, and 
(c) <j> = — 2 for the plus sign in eq. (2.16). The left figures correspond to past solutions while right figures 



correspond to future solutions. Each region corresponds to (I) nonsingular solutions with D — > 2, (II) 
singular solutions with determinant singularity (D = 0), (III) singular solutions with D — » — oo, (IV) 
nonsingular solutions with D — > +oo, respectively. The black color indicates prohibited regions in the 
initial condition space. 
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FIG. 3: The minimal initial values of d> which allow solutions nonsin- 



gular in future in the dilaton-driven case for the plus sign in eq. (2.16) 



B. Minus sign of eq. (2.16) 



In the previous subsection solutions nonsingular in future meet the determinant singularity in past, due to the fact 
that the determinant evolves from the asymptotic value D — 2 toward the determinant singularity. For the minus 



sign of eq. (2.16), the determinant can increase toward the past, implying that the trajectories may not cross the 
determinant singularity. In fact we will show that the past singularity is not necessarily of the determinant-type. 

Let us first analyze the asymptotic behavior of the determinant when all expansion rates diverge as eq. ( 3.1 C| ) . In 
eq. (2.16) when the 2(pq + qr + rp) term is positive and much larger than (2Apqrf) 2 , one has (f> ~ — \/2(pq + qr + rp) ~ 
—y6po/(t — t s ). Therefore the dilaton grows as cf) ~ — V6po\n (t — t s ) for t — t s — > +0, in which case the asymptotic 
form of the determinant yields 



D 



2 + 6V6X P 2 {t - £ s ) 2 (^po-i) + 39A 2 p 4(t _ ts )4(V6 P0 -i) 



18V6\ 3 p 6 ( 



(3.16) 



This indicates that the determinant grows infinitely (D — > +oo) for po < 1/ V6, while it decreases toward a finite value 
D = 2 for po > 1/ y6. The latter case corresponds to the one where the average expansion rates are large initially. 

When the condition (2Apqrf) 2 ^> \2(pq + qr + rp)\ is satisfied, the asymptotic form of <f> is given by <j> ~ 
— 6\pQe~ 2r ^ / (t — t s ) 3 , which yields e 2 ^ ~ 6Apg/(i — ts) 2 - Then one has the following positive finite value of the 
determinant in the limit of t — t s — > +0: 

13 7 5 



D~2 + 



(3.17) 



I2p 2 36p4 432^- 

Whether the asymptotic form is given by (3.16) or ( |3.17j ) depends on the initial values of the expansion rates and the 
dilaton. In the latter case <j> is generally small so that the condition (24pqrf) 2 ^> 2(pq + qr + rp) is fulfilled. In both 
cases we can expect that the trajectories do not cross the determinant singularity in past. 

We show in Fig. || the density plot of the asymptotic behavior of the determinant for h = 0.16 at t = 0. When 
4> = at t = all trajectories satisfying the constraint (2.15) are nonsingular in future with D = 2. In this case the 
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past solutions are dominated by the singularity with asymptotic positive determinant as found from the left panel of 
Fig. H (a). As one example we plot in Fig. [| (a) the evolution of p, q, r, <fi, and D for <fi — at t = with anisotropy 
parameters a — 0.05 and (3 — 0.05. This belongs to the class ( 3.1 6| ) with a. < 1/V6 where all expansion rates diverge 
with determinant D — > +oo. The determinant is always positive and continues to grow pastwards. The left panel 
of Fig. U (a) shows that there exist some past solutions which meet the determinant singularity when anisotropies 
parameters are large. In this case the determinant grows until some moment of time pastwards, after which it begins 
to decrease toward the determinant singularity. This behavior can be understood that large anisotropies prevent all 



expansion rates from evolving almost similarly as eq. (3.10). However this region is typically small for the positive 
values of 6 at t = 0. 



With the decrease of the initial cj), the constraint ( 2.15 ) gives wider allowed parameter ranges in (a, /3)-planc. Let 
us consider the case cf) = — 1 and h — 0.16 at t = 0. As is found from Fig. || (b) we have additional regions with 
the D — — oo singularity in both past and future. Nature of singularities around the isotropic point is similar to the 
(f> = case explained above. For the solutions nonsingular in future, the past singularity is either the type of D > 
or D — 0. Fig. H (b) is the latter case where the trajectory comes regularly from the asymptotic future (D = 2) and 
meets the determinant singularity in past. This evolution is similar to the case of large anisotropy parameters in Fig. |5| 
(a). From Fig. ^ (b) the solutions singular in future (D — — oo) correspond to, in past, the singularity with D = — oo 
or the determinant singularity. The difference between two cases is whether the determinant is always negative [like 
Fig. [l] (b)] or it passes through D = 0. 

The situation becomes somewhat different with <j> being decreased further. Fig. || (c) indicates that nonsingular 
trajectories with D — > +oo appear in future around the isotropic point for (f> < — 2 at t = 0. This is the case all 
expansion rates are finite as described by eq. ( |3.5| ). However these solutions are singular in past with all expansion 
rates being infinite. In the region V shown in Fig. || (c), the condition (2Apqrf / ) 2 > \2(pq + qr + rp)\ is typically 
satisfied at t = 0. Therefore the determinant tends to approach the finite value ( |3.17 ) in past. We show in Fig. [| (c) 
the evolution of the system for a — —0.02, (3 = 0.02, cj) = —2, and h = 0.16. While the evolution of the expansion rates 
in Fig. [| (c) looks similar as in Fig. ^ (a), the behavior of the determinant is different in both past and future. In Fig. ^ 
(c) the determinant decreases from infinity to a finite positive value toward the past. We have numerically found that 
two terms (24pqrf) 2 and 2{pq + qr + rp) become comparable during the evolution. Therefore the asymptotic value 
of D is somewhat different from eq. ( 3.17 ), but it is still a finite positive value. We find from Fig. ||| (c) that future 
nonsingular solutions with either D = 2 or D = +oo correspond to the singularity ( |3 .10 ) with positive determinant 
in past. Figs. || (b) and (c) also suggest that as (j> decreases the ordinary singularity with D — > — oo appears in both 
past and future for large anisotropy parameters, 

Compared to the case of the previous subsection, the range of nonsingular solutions in future is not so narrow due 
to the presence of nonsingular trajectories with D — > +oo. However this solution is not appropriate to lead to our 
present Universe due to the dominance of the quadratic curvature term. In addition these solutions typically approach 
the ordinary singularity (3.10) with all expansion rates infinite in past. This property is different from the plus sign of 
eq. (2.16) where the past singularity for solutions nonsingular in future corresponds to the determinant-type. In both 
cases we have found that nonsingular cosmological solutions in both past and future do not exist for the dilaton-driven 
case for wide ranges of the parameter space (10~ 5 < h < 1 and — 15 < <j> < 15). However the situation is changed in 
the modulus-driven case as we will analyze in the next section. 
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(b) 





FIG. 4: The evolution of p, q, r, <f), and D in the dilaton-driven case for the minus sign in eq. ( 2.16 ). 
Each figure corresponds to initial conditions with h — 0.16 and (a) <f> = 0, a = 0.05, j3 = 0.05, 
(b) (j> = -1, a = 0.15, 13 = 0.05, (c) </> = -2, a = -0.02, /3 = 0.02, respectively. 
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FIG. 5: Nature of singularities in (a, /3)-plane with initial conditions, (a) (f> = 0, (b) cf> — 
— 1, and (c) <f) — —2 for the minus sign in eq. ( 2.16 ). The left figures correspond to past 
solutions while right figures correspond to future solutions. Each region corresponds to 
(I) nonsingular solutions with D — > 2, (II) singular solutions with determinant singularity 
(D — 0), (III) singular solutions with D — > — oo, (IV) nonsingular solutions with D — > oo, 
(V) singular solutions with positive determinant, respectively. The black color indicates 
prohibited regions in the initial condition space. 
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IV. MODULUS-DRIVEN CASE 



In the modulus-driven case the function /(</>) in eq. ( |2.1| ) is expressed as P jK| , |l2| -]l-j 

1 



16 



(4.1) 



where the coefficient S is determined by the 4-dimensional trace anomaly of the N = 2 sector. Here the function £(</>) 
is defined by 



£((/>) = In [2eV(^)] 



(4.2) 



where r]{ie^) is the Dedekind 77 function. Since £(0) is well approximated as £(0) ~ — (27r/3)cosh</> |14[ ], the function 
/ takes the form 

n5 , jl . 



/(£> 



48 



(4.3) 



When is largely negative (|(/>| S> 1), eq. ( |4.3| ) reduces to the form (3.1) by setting <$ = 3A/7T. Therefore when 5 > 
solutions nonsingular in both past and future do not exist in this case. However nonsingular cosmological solutions 
have been found for negative values of S Hereafter we shall focus on the negative S case (setting 5 = —1 for 
simplicity) using the function /(</>) given by (|4.3|). 

When the Gauss-Bonnet term dominates in eqs. (2.7)- (|2.10 ), the asymptotic solution is similarly given as eq. (pi]), 
i.e., 



P = dm 



q = C2\t\ 



C3\tl 



6 ±51n|i|, C1C2C3 = sign(t)- 



AttS 



(4.4) 



The past asymptotic solutions correspond to C1C2C3 > for negative S. This sign is different from the dilaton-driven 
case with positive A which we already analyzed in the previous section. Note that the determinant is divergent 
(D -> +00) for \t\ -> 00. 

For another asymptotic solution where the effect of the Gauss-Bonnet term is negligible, the evolution of the 
background is given by 



Ci+C 2 + c 3 



Q = ca|t| 
= sign(t), 



r = c 3 \t\ 1 , = 0o + c 4 ln|t|, 



2 , 2 1 2 
Cj + c 2 + c 3 



2(cic 2 + c 2 c 3 -I- C3C1) . 



(4.5) 



Therefore one has /' oc |t|' C4 ' , cj> oc |t| _1 , and p,q,r cx lil" 1 with | C4 1 < 1 in eq. (2.14), in which case the solutions are 
nonsingular with an asymptotic value of the determinant, D = 2. 

For the nonsingular cosmological solutions found in ref. |J , does not change its sign jl(J . Due to the symmetric 
structure of the function (4.3) with respect to <f> = 0, we will consider the positive <f> case where the modulus continues 
to grow from past to future. For negative 4> the analysis is essentially the same by changing (f> to —<j). When the 
solutions are singular, they meet the determinant singularity [see eqs. (|3.6[) and ( |3.7|) 1 or the ordinary singularity [see 
eqs. (3.8)-( 3~10| )l. Let us consider the asymptotic behavior of the determinant for the ordinary singularity. In the 
cases ( 3.8 )-( 3.10 ), possible asymptotic behavior of the determinant can be summarized as 



Case(i) 

Case(ii) 
Case(iii) 



D 

D 
D 



3(i-i s ) 2 

3(i-i. 5 ) 4 
0, 4> ~ - 



( q °- r i)' 



with 



1 go + r 



const . 



-00, 



with 



const . 



t-t. 



(4.6) 

(4.7) 
(4.8) 
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where f' s = (7r<5/48)(e* s — e~^ s ) with <f> s being a constant. Although the case (ii) is rare, we have numerically checked 
that this asymptotic solution certainly exists for the modulus-driven case. The asymptotic (iii) corresponds to the 
isotropic case, as is similar to the dilaton-driven cosmology. It is easy to verify that this trajectory is possible for 
positive 4> but impossible for negative (j>. In what follows we shall analyze how the evolution of the expansion rates 
and the determinant are different from that of the dilaton-driven case paying particular attention for nonsingular 
trajectories. 

Nonsingular asymptotic solutions for (f> > can be described as 4> = 4>q — 51n|t| in past [see eq. (fO|)] and 4> — 
4>q + C4 In I i I with < C4 < 1 in future [see eq. (4.5)]. Therefore <j) continues to grow from asymptotic past starting 
from large negative values of (j) toward the future. Such examples are plotted in Fig. || (a) and (b). In these cases the 
determinant evolves from D = +00 (past) to D — 2 (future) without passing through the determinant singularity. In 
the dilaton-driven case when the past trajectories are nonsingular with D = +00 they inevitably meet the determinant 
singularity in future [see Fig. || (b) and (c)]. This is mainly due to the fact that /' is negative in the dilaton case 
while its sign is different in the modulus case for 8 < and <fi < 0. Therefore the determinant ( 2.14 ) is dominated by 
positive terms in the modulus case, which provides a way not to pass through the determinant singularity. Namely 
negative d is crucial for the existence of nonsingular solutions. 

In Fig. [7] we show past and future asymptotic properties in three different cases. Note that we have defined a = a/h 
and [3 = (3/h in order to compare the cases where the average expansion rate h is changed. When h = 0.05 and 
<p = —5, nonsingular solutions in past (D = +00) are not singular in future with determinant, D — > 2 [see Fig. ^ (a)]. 
The Universe exhibits superinflation with growing expansion rates until the graceful exit around t = [see Fig. ^ (a)]. 
Notice that we have pqr > in asymptotic past, as predicted by eq. (4.4) for negative S. The expansion rates begin 
to decrease after the graceful exit, whose asymptotic solutions in future are given by eq. Q. In Fig. | (a) we find 
that the future solution corresponds to the Kasner-type where the Universe is contracting in one direction. 

Fig. ^ (a) indicates that some trajectories which are nonsingular in future cross the determinant singularity in past. 
The evolution of the background is similar to Fig. [j] (a) which we already analyzed in the dilaton case. We also find 
from Fig. ^ (a) that when anisotropy parameters are large the solutions meet the ordinary singularity with D = — 00 
in both past and future. This is the case of (3.8) or (3.9) where at least one expansion rate diverges as plotted in 
Fig. [I] (b). Although in Fig. [7] past nonsingular solutions (D — > +00) do not meet the determinant singularity in 
future, we have checked that this singular behavior occurs for smaller values of <\> as shown in Fig. [l] (c). These results 
imply that the property of singularities is similar to the dilaton-driven case described in Sec. Ill A. 

For larger initial values of cj>, the allowed region can be approximately described as a 2 + 1 < 1. When h = 0.05 
and (j> = 2 shown in Fig. [7] (b), all future asymptotic solutions are nonsingular with determinant, D — > 2. However we 
find from the left panel of Fig. (b) that the region of past nonsingular solutions gets smaller relative to Fig. |?] (a). 
For large anisotropy parameters, the determinant tends to decrease pastwards from asymptotic future with D = 2, 
thereby resulting in the singularity with D = 0. This is not the case for the nonsingular trajectories with small 
anisotropy parameters. One example is plotted in Fig. ^ (b). Although the determinant decreases from t = +00 to 
t ~ —6, it begins to grow toward the asymptotic past before crossing the determinant singularity. From Fig. ^ (b) 
we find that the trajectory comes regularly from asymptotic past and connects to the Friedmann-type branch where 
the Universe is expanding in all directions. When anisotropy parameters are small and belong to the region ( |2.6| ), the 
future solution is of the Friedmann-type. 
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FIG. 6: The evolution of p, q, r, <f>, and D in the modulus-driven case. Each figure corresponds to 
initial conditions with (a) <j> = -5, h = 0.05, a = a/h = -0.5, (3 = /3/h = 0.5, (b) = 2, h = 0.05, 
& = 0.3, = -0.2, (c)(f> = 2,h = 0.1, & = 0.5, (3 = -0.5, respectively. 
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FIG. 7: Nature of nonsingular and singular solutions in (<5,/3)-plane with initial 
conditions, (a) h = 0.05, (f> = -5, (b) h = 0.05, <j> = 2, and (c) h = 0.1, = 2. 
The left figures correspond to past solutions while right figures correspond to future 
solutions. Each regions correspond to (I) nonsingular solutions with D — > 2, (II) 
singular solutions with determinant singularity (D = 0), (III) singular solutions with 
D — » — oo, (IV) nonsingular solutions with D — > oo, respectively. The black color 
indicates prohibited regions in the initial condition space. 
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FIG. 8: The regions of nonsingular and singular solutions for the 
modulus-driven case in the initial condition space of h and 4>. 



When the average expansion rate h gets larger, the region of nonsingular cosmological solutions becomes smaller. 
We plot in Fig. fj] (c) the density plot of nature of the determinant for h = 0.1 and <j> = 2. The past solutions are 
dominated by the determinant singularity, whereas the future solutions are nonsingular as is similar to Fig. [t] (b) . In 
this case nonsingular trajectories in both past and future exist in only small parameter ranges around the isotropic 
point. 

It is also worth mentioning that the ordinary singularity (3.10) with asymptotic determinant ( |3.16|) and (3.17) for 
the minus sign of eq. ( 2.1 6| ) does not appear in the present case. When <f> increases toward the past, |/'| also grows 
in the modulus case. Therefore the condition (2Apqrf) 2 <C \2{pq + qr + rp)\ is not satisfied, implying that case of 
(3.16) does not occur. When (24pqrf) 2 3> \2(pq + qr + rp)\, there are two possibilities for the asymptotic form of 0, 
one of which is cb ~ ir8(e^ — e~ ( >')p^ s /{t — t s ) 3 and another is <f> ~ C[t — t s ) where C is a constant. In the former case 
it is easy to show that asymptotic solutions do not exist by integrating <fi ~ n8(e^ — e _ ^)po/(i — t s ) 3 [Note that this 
is possible in the dilaton-driven case with determinant (3.17)]. In the latter case the determinant approaches D = 
as in eq. ( ]4.S| ), 

We have done numerical simulations for other cases varying the values of h and <p. When h is large, <j) is required 
to be small for the existence of nonsingular cosmological solutions. This property is found in Fig. ^| where we plot 
the regions of nonsingular and singular solutions in the (h, 0)-plane. Nonsingular trajectories come regularly from 
asymptotic past with D = +oo where the quadratic curvature term is dominant [see eq. (4.4)], and smoothly connect 
another nonsingular branch with D = 2 where the Gauss-Bonnet term is negligible [see eq. (4.5)1. F° r singular 
solutions nature of singularities is found to be similar to that of the dilaton-driven case discussed in Sec. Ill A. 



V. SUMMARY AND DISCUSSIONS 



In this paper we have analyzed past and future asymptotic regimes in Bianchi type-I string-inspired cosmological 
models in the presence of a Gauss-Bonnet curvature invariant. If the loop correction is not taken into account, one has 
no-go results that the initial big-bang singularity can not be avoided. The Gauss-Bonnet term allows the existence of 
nonsingular cosmological solutions, depending on the theories we adopt. We investigated two gravity theories, viz., 
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dilaton- and modulus-driven cosmologies. In the former case the dynamics appears to depend significantly on the 
sign of rhs of cq. (|t|). Hence, we treated possible signs in eq. ( |2.16|) separately and constructed three pictures of 
cosmological evolution: 



(a) Dilaton-driven cosmology with plus sign in eq. (2.16) 

(b) Dilaton-driven cosmology with minus sign in eq. ( [2.16 ) 

(c) Modulus-driven cosmology with S < . 

As the quadratic curvature corrections may provide violations of strong and week energy conditions Jl3| , p^[ , the 
Bianchi type-I Universe can recollapse in high-curvature regime. This fact gives us a variety of possible types of 
trajectories with different past and future asymptotics. Some trajectories can not leave high-curvature regime, some 
reaches the low-curvature future attractor. 

In the dilaton-driven cosmology the 4-dimensional string coupling, A, is required to be positive. This forbiddens 
the existence of nonsingular cosmological solutions even when the Gauss-Bonnet term is taken into account. In the 
case (a) nonsingular trajectories in future correspond to the low-curvature solutions with determinant D — > 2 where 
the Gauss-Bonnet term is negligible. These trajectories meet the singularity where the determinant of the system 



vanishes at finite past. At this determinant singularity, p,q,f,4> diverge in eq. (2.12) while p,q,r,<j) stay finite. This 



kind of singularities restricts the presence of nonsingular solutions in an anisotropic background. There exist other 



kind of singularities (we call ordinary singularities) where at least one expansion rate diverges [see eqs. (3.8)-(3.10)] 



When the singularity (3^8) or (3J3) appears at finite past or future, the determinant approaches D = — oo. We also 
find some trajectories which are nonsingular in asymptotic past with determinant, D — > +oo. These solutions appear 
for large anisotropy parameters, whose existence, to our knowledge, was not discovered previously. In the case (a) 
these solutions are found to meet the determinant singularity in future (see Fig. 0). In the case (b) there appears 



the ordinary singularity of eq. (3.10) where all expansion rates diverge. In this case the determinant is divergent 



or approaches a positive constant value, depending on the initial conditions of 4> and h [see ( |3.16j ) and (3.17)]. The 



solutions nonsingular in future typically correspond to this type of singularity or the determinant singularity in past 
as found by Fig. [g. 



In the modulus-driven case the coupling, 5, can take either positive or negative value in eq. (4.3). When <5 > 
solutions nonsingular in both past and future ar not found as is similar to the dilaton-driven case. For negative values 
of 5, however, there exist nonsingular cosmological solutions where two branches of superinflation and decreasing 
curvature can be joined to each other. In this case the solutions come regularly from asymptotic past (D = +oo) 
with determinant being decreased. The determinant approaches the future asymptotic value D = 2 without crossing 
the determinant singularity. This is the main difference from the dilaton-driven case where the past nonsingular 
trajectories inevitably meet the determinant singularity in future. When the solutions are singular we find that 
nature of singularities is similar to the dilaton-driven case discussed in Sec. Ill A. 

Our numerical investigations show that not all combinations of past and future asymptotics are possible. We pay 
particular attention to nonsingular past and low-curvature nonsingular future regimes (the latter should correspond 
to our present Universe). For these regimes we get the following results: 

• For the dilaton-driven case trajectories with nonsingular past asymptotic meet the determinant singularity in 
future. 

• For cases (a) and (c) trajectories with low-curvature future asymptotic meet the determinant singularity in past, 
when the past solutions are singular. 
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The first property tells us that the negative sign of the coupling constant is essentially important for constructing a 
purely nonsingular string cosmological models. This result is known for FRW Universe [p|JlO|l and, hence, is still valid 
in anisotropic Bianchi I case, despite the existence of a past nonsingular regime in the dilaton-driven cosmology. The 
second property indicates that the determinant singularity may play a crucial role when we try to trace back in time 
evolution of the our present Universe in anisotropic background. 

Recently string-inspired cosmological models which can avoid the big bang singularity have received much attention 
p0| pl| together with the proposal of the Ekpyrotic universe |2j|. In those cases the quantum loop corrections or 
the higher-order derivatives play important roles to determine the dynamics before the graceful exit. It is certainly 
of interest to extend our analysis to more complicated models such as the multi-field case in the presence of the 
Gauss-Bonnet term. In fact while there exist nonsingular solutions in the single-held modulus-driven case considered 
in this work, density perturbations generated by the fluctuation of modulus exhibits blue-spectra with a spectral tilt 
n = 10/3 |2(|g7]| . This contradicts with the observational supported flat spectra with n ~ 1. In the multi-field case, 
however, it may be possible to produce almost scale-invariant spectra if a light scalar field such as axion generates a flat 
isocurvature perturbations during superinflation |2^-|30[| or if the axion is nonminimally coupled to the dilaton with 
some potential [ [3l| . We leave to future work to construct such nonsingular cosmological models which are consistent 
with observations. 
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